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Abstract 

We apply a quark interchange model to spin-dependent and exotic 
meson-meson scattering. The model includes the complete set of stan- 
dard quark model forces, including OGE spin-orbit and tensor and scalar 
confinement spin-orbit. Scattering amplitudes derived assuming SHO and 
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Coulomb plus linear plus hyperfine meson wavefunctions are compared. In 
1=2 TTTT we find approximate agreement with the S-wave phase shift from 
threshold to 1.5 GeV, where we predict an extremum that is supported by 
the data. Near threshold we find rapid energy dependence that may rec- 
oncile theoretical estimates of small scattering lengths with experimental 
indications of larger ones based on extrapolation of measurements at mod- 
erate fcyr^. In PsV scattering we find that the quark-quark L-S and T forces 
map into L-S and T meson-meson interactions, and the P-wave L-S force is 
large. Finally we consider scattering in J^*^"-exotic channels, and note that 
some of the "Deck effect" mechanisms suggested as possible nonresonant 
origins of the 7ri(1400) signal are not viable in this model. 



I. INTRODUCTION 

The determination of scattering amplitudes between pairs of mesons is an interesting 
problem in strong QCD. It is also a complicated problem, because both qq annihilation to 
s-channel resonances and "nonresonant" scattering are important effects, and it is often 
difficult to separate the various contributions. However by specializing to annihilation- 
free channels such as 1—2 nn and 7rp, 1=3/2 Ktt, KN and NN, one may study nonresonant 
scattering in relative isolation. The determination of resonance parameters, reaction 
mechanisms, and many other aspects of hadron physics are complicated by the presence 
of nonresonant scattering, which is treated as an (often poorly understood) initial-state 
and final-state rescattering effect. Developing an accurate description of nonresonant 
scattering would help clarify many other aspects of hadron physics. 

A further interesting possibility is that sufficiently attractive nonresonant scattering 
may lead to weakly bound hadron-hadron or multihadron states, as does happen in nuclei 
and hypernuclei. We may also find a rich spectrum of meson-meson bound states, the 
study of which will extend nuclear physics into the largely unexplored field of "mesonic 
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nuclei" or "molecules" 1T|-§|. 

An understanding of PsPs, PsV and other meson-meson scattering amplitudes is also 
important for the interpretation of non-QCD processes such as nonleptonic weak decays, 
since these show evidence of important hadronic final state interactions. The Al=l/2 rule 
is a well known example. Similarly, a recent study of D and decays to KKvr found 
that the Dalitz plots are dominated by two-meson isobars, including 07r, K*K + h.c. and 
Kq(1430)K + /i.e., and complex relative amplitudes are required to describe the D"^ Dalitz 
plot. Without final state interactions one would expect relatively real couplings to these 
final states. 

One finds a surprising variety of approaches to strong hadron-hadron scattering in the 
literature. There are many studies using effective hadronic lagrangians, such as the "chiral 
perturbation theory" description of the PsPs sector. Although this method is convenient 
because it uses perturbative QFT techniques, it is incomplete in that it takes effective 
lagrangian vertex strengths from the data; one should be able to calculate these hadronic 
couplings directly from quark-gluon forces. 

Second, there are studies that model the low energy hadron-hadron scattering mech- 
anism, which include the apparently dissimilar meson exchange and quark-gluon descrip- 
tions of hadronic forces. Meson exchange models are again attractive for their simplicity, 
since they use perturbative QFT techniques to determine scattering amphtudes. This 
approach has been elaborated in greatest detail in models of the NN force 0, in which 
a large number of meson exchanges is assumed. With this large parameter space a good 
description of this interaction is possible, although there is a concern that one may be 
parametrizing other scattering mechanisms in addition to t-channel meson exchange. Al- 
ternatively, one may calculate hadron-hadron forces directly from the fundamental quark- 
gluon interaction, using quark model hadron wavef unctions. This approach has also seen 
its most detailed development in studies of the NN interaction 0, and is most success- 
ful in describing the short-ranged repulsive core. Maltman and Isgur [^] also found a 
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physically reasonable intermediate ranged attraction from a color van der Waals effect in 
the quark-gluon approach, which is not equivalent to the usual tttt or a meson exchange 
explanation of this force. The quark description of hadron-hadron interactions is com- 
plicated by the combinatorics of matrix elements between quark bound states, but has 
the advantage that it can easily be extended to a wide range of spin and flavor channels 
through a simple change of the external hadron wavefunctions. 

A third promising approach is to infer hadron scattering amplitudes from LGT. To 
date LGT has seen little application to scattering problems because of the difficulty of 
treating systems that are not in their ground states. Estimates of the 1=0 and 1=2 tttt 
scattering lengths have been obtained by exploiting a theoretical relation to finite-size 
effects 0, and more recently very interesting results for nuclear physics potentials in the 
BE system were reported P]. In future it may be possible to improve hadron scattering 
models through comparisons with similar "LGT data" . 

In this paper we are concerned with the derivation of meson-meson scattering ampli- 
tudes from quark-gluon forces. We derive meson-meson scattering amplitudes at lowest 
order in the quark-quark interaction, which leads to a quark interchange model described 
by "quark Born diagrams" P,pi]|. Since the quark-quark interaction is already well es- 
tablished from hadron spectroscopy, our predictions have little parameter freedom. In 
previous work we and others (usually assuming OGE hyperfine dominance) have shown 
that this approach gives a reasonably accurate description of S-wave scattering in a wide 
range of channels without qq annihilation, including 1=2 tttt Q, 1=3/2 Ktt [|rT|, 1=0,1 
KN PI, 1=0,1 BB (compared to LGT data), and the NN repulsive core [g. This 
approach has also been applied to TtJ/ip [ll4| , p!5| and other reactions relevant to heavy ion 
collisions, where the experimental low energy cross sections are as yet unclear. 

The principal new contribution of this paper is a detailed analytical derivation of the 
meson-meson scattering amplitudes that follow from the complete quark-quark interac- 
tion, including color Coulomb, linear scalar confinement, OGE spin-spin, OGE spin-orbit. 
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OGE tensor and linear spin-orbit forces. As a future application of these results, one might 
hope to clarify the relationship between meson exchange and quark interchange models 
by a detailed comparison of the spin dependence of hadron-hadron scattering amplitudes, 
which we expect to be sensitive to the details of the scattering mechanism. 

Here we consider both pseudoscalar-pseudoscalar (PsPs) and pseudoscalar-vector 
(PsV) scattering. The former is a "standard benchmark" for meson scattering mod- 
els, because 1=2 tttt low energy scattering has no s-channel resonances and has been the 
subject of many experimental phase shift analyses. Although we find reasonable agree- 
ment with S-wave 1=2 tttt scattering, this channel has no spin degree of freedom, and so 
cannot be used to test the characteristic spin dependences predicted by the quark model's 
OGE and linear scalar confinement forces. 

We find in contrast that PsV is an excellent theoretical laboratory for the study of 
spin-dependent forces, as it can accommodate both meson-meson spin-orbit and tensor 
interactions. The spin-dependent forces at the meson-meson level are closely related to 
the corresponding terms in the quark-quark interaction in our approach. Although the 
study of PsV scattering is essentially a theoretical exercise at present, these phase shifts 
are accessible experimentally, for example through measurement of the relative S and D 
final state phases in bi txuj. Thus it should be possible to measure PsV phase shifts 
from resonance decays to multiamplitude PsV final states. 

Before we proceed to our detailed results, we note that some work has already appeared 
on meson-meson scattering in PsV systems. Numerical results for many light S-wave PsV 



meson channels were previously reported by Swanson [|T0[ using a similar quark model 
approach that incorporated OGE spin-spin and linear confinement forces. Theoretical 
results for PsV scattering (vrp in particular) in a meson exchange model were published 



by Janssen et al. |T6| and Bockmann et al. |T^, assuming vr, vector and ai exchange. Since 
the pTTTT, aipvr and ptuvr vertex strengths are relatively well established, it was possible to 
evaluate these scattering amplitudes numerically. These papers did not consider the exotic 
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1=2 channel, so a direct comparison with our quark model PsV results is not possible at 
present. 

II. MESON-MESON T-MATRIX 

A. General T-matrix formula 

We approximate the full hadron-hadron scattering amplitude by a single (Born-order) 
matrix element of the quark-quark interaction Hamiltonian Hi. Since Hi is T°'T'^ in color, 
one must then have quark line rearrangement to have a nonvanishing overlap with two 
color-singlet mesons in the final state. In {qq) — (qq) scattering there are four independent 
Born-order diagrams, which we label according to which pair of constituents interacted; 
these are "transferi" (Tl), "transfer2" (T2), "capturei" (CI) and "capture2" (C2), which 
are shown in Fig.l. In the special case of identical quarks and identical antiquarks, which 
is relevant here, there is a second set of four "symmetrizing" diagrams Tlsymm ■ ■ ■ C2symm, 
which are identical to Tl. . .C2 except that the quark lines are interchanged rather than 
the antiquark lines. 

The hadron-hadron T-matrix element Ty, for each diagram can conveniently be written 
as an overlap integral of the meson wavefunctions times the underlying quark Tfi. These 
overlap integrals (specializing Ref. [|T^] to the case of equal quark and antiquark masses) 
are 

T^J^\aB ^CD) = 

JJ SqSp ^*c{2p + q-C) ^*j^{2p-q-2A-C) 

Tf,{q,p,p-A-C) ^Ai2p-q-A) ^B{2p + q-A-2C) , (1) 

T^r\AB CD) = 

JJ SqSp ^*c{-2p + q + 2A-C) ^*o{-2p-q-C) 
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(2) 



T}^'>{AB CD) = 

JJ d^qSp ^};{2p + q-C) $^(2p-g-2l-(7) 

TfMP.-P + C) ^Ai2p-q-A) <^Bi2p-q-A-2C) , (3) 

Tjf^\AB CD) = 

JJ SqSp ^*c{-2p + q + 2A-C) $^(-2p-g-C) 

TfMP.-p-C) ^Ai-2p + q + A) ^s(-2p + q + A-2C) . (4) 

The quark Tfi has momentum arguments Tfi{q,pi,p2), which are defined in Fig. 2. In 
this paper we will evaluate these overlap integrals with standard Gaussian quark model 
wavefunctions and the quark Tfi for the complete set of OGE color Coulomb, linear 
scalar, OGE spin-spin, OGE and linear scalar confinement spin-orbit and OGE tensor 
interactions. These interactions are given in App.A. 

B. PsPs Scattering 

1. 1=2 vrvr T-matrix 

We specialize the general problem of PsPs scattering without qq annihilation to 1=2 
TTTT because many experiments have published phase shift analyses of this channel. The 
other TTTT channels have large s-channel qq annihilation contributions. The full 1=2 vrvr 
Born-order T-matrix element is determined by adding the individual contributions of 
App.B, with PsPs spin matrix elements given in App.C, part C2. There are also flavor 
and color factors for each diagram and an overall "signature" phase of (—1), and a second 
set of "symmetrizing" diagrams for identical quarks and identical antiquarks, as discussed 
in detail in Ref. [^. On summing these contributions we find 
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where ia,c{x) is an abbreviation for the confluent hypergeometric (Kummer) function 
iFi(a;c; x). 

The three separate expressions above are the OGE spin-spin, color Coulomb and linear 
confinement contributions respectively. The Q± terms come from the transfer diagrams, 
and the remaining, isotropic, terms come from the capture diagrams. The spin matrix 
elements of the spin-orbit and tensor terms vanish identically in the PsPs channel. 

Since Q± = C ± A and \A\ = \C\, one can equivalently write this amplitude as a 
function of the CM momentum and scattering angle using = 2A^{1 ± /i), where 
H = cos{6ac)- The Bose symmetry required for this vrvr scattering amplitude is evident. 



2. 1=2 -K-K Phase Shifts 

We may derive the elastic Born-order 1=2 tttt phase shifts from Eq.(|^), using the 
relation between phase shifts and the T-matrix given in App.D, especially Eq.(D17), and 
the integrals in App.G. The result we find for the S-wave is 



r7=2 TTTT 

■'o 



32 X 



-2x 



37/2 



-4x/3 



OGE spin-spin 

kE^f ( - ^i(fi,i (-2a;) - e-^-) - fj^ fij (-2x/3) e'^^/^) OGE color Cou. 
kE^j^ (^^l{{^ i{-2x) - + ^ f2j(-2x/3) e'^^/^^ . linear conft. 

(6) 

-> 2 

where we have introduced x = A /4/?^. The total Born-order S-wave phase shift is the 
sum of these three contributions. 

This S-wave phase shift is shown in Fig. 3 with our standard quark model parameter set 
a, = 0.6, p = 0.4 GeV, m = 0.33 GeV and b = 0.18 GeV^. We also use = 0.138 GeV 
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throughout. This confirms that the color Coulomb and linear confinement interactions 
make relatively small contributions to the 1=2 tttt S-wave at moderate energies. The 
weakly repulsive linear confining interaction in 1=2 tttt near threshold was previously 



found numerically by Swanson ||T^ . 

One might be concerned about the approximation of using SHO wavefunctions, es- 
pecially at higher energy scales where there should be strong short-distance components 
in the pion wavefunction due to the attractive spin-spin hyperfine interaction. To test 
the sensitivity to SHO wavefunctions we evaluated the 1=2 tttt scattering amplitudes 
and phase shifts numerically using Coulomb plus linear plus hyperfine qq wavefunctions 
and Monte Carlo integration of the real space integrals corresponding to the T-matrix 
integrals (^^. As usual this requires a "smearing" of the contact hyperfine term, 
6{x) e^'^^'^^/vr^/^a^/^, to allow solution of the Schrodinger equation with an attrac- 
tive delta-function interaction. In the literature the inverse smearing length is typically 
taken to be cr ~ 1 GeV. (A calculation of 1=2 tttt scattering with this interaction and 
a = 0.7 GeV was reported previously by Swanson [jl^].) With our standard light-quark 
parameter set = 0.6, b = 0.18 GeV^ and m = 0.33 GeV, we found that fitting the 
Mp — splitting required a value of a = 0.86 GeV. To illustrate the dependence of the 
scattering amplitude on this parameter, in Fig. 4 we show the 1=2 tttt S-wave that follows 
from our standard quark model set (a^, b, m), with a = 0.7, 0.8 and 0.9 GeV. Clearly the 
predicted phase shift is rather similar to the SHO result of Fig. 3, although the effect of 
short- distance peaking in the vr wavefunction is evident above M^^ ^ 1 GeV. 

We also show most of the higher statistics experimental results for the 1=2 tttt S-wave 
phase shift in Figs. 3 and 4. The references shown are Colton et al. Durusoy et al. 
||T9| (OPE extrapolation, solid; OPE + DP form factor, open, slightly displaced in x for 
visibiity), Hoogland et al. |^ (extrapolation B), and Losty et al. [|2T|. Prukop et al. 
found a wide range of results from three different off-shell extrapolations, so we simply 
quote their fitted scattering length below. 
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Clearly there is already reasonable agreement with the experimental S-wave phase 
shift at lower energies without fitting the quark model parameters. The model predicts a 
rather dramatic extremum in this phase shift near M^^ = 1.5 GeV, which is unfortunately 
beyond the limiting invariant mass of most of the experimental studies. There are some 
measurements of this phase shift at higher invariant mass with lower accuracy due to 
Durusoy et al. ||19|, which are also shown in the figure. The Durusoy et al. results 
support our predicted extremum near Mt,t, = 1.5 GeV; indeed, their phase shift above 
Mtttt = 1.5 GeV appears to fall even more rapidly than we predict. 

We have investigated optimal parameter fits of the S-wave phase shift formula Eq. (6) 
to the data, but we find that these are rather unstable because the color Coulomb and 
linear confinement contributions are small and are qualitatively similar functions. In any 
case the Durusoy et al. data and Fig. 4 show that the hyperfine smearing distance a is 
an important parameter, and this will not be well determined until accurate phase shift 
data becomes available at higher invariant mass. An accurate measurement of 1=2 tttt 
scattering amplitudes near and above M^^^^ = 1.5 GeV would clearly be very useful as a 
test of this and other models of meson-meson scattering. 



3. 1=2 TTTT Scattering Lengths 



The 1=2 TTTT scattering length is defined by Oq ^ = limfc^^o 
find from Eq.(6) are 



1 = 2 TTTT 



/kj,. The results we 



.7=2 



_2_ 

'32 



2» 



< ^2 I 5172 1 



2 



M 



S-S 
Cou. 
lin. 



(7) 



and their numerical values with our standard quark model parameters set are 



10 



J=2 



' -0.085 [fm] S ■ S 
-0.007 [fm] Cou. 

(8) 

-0.017 [fm] lin. 
^ -0.109 [fm] total. 

The Coulomb and linear contributions were independently checked by Monte Carlo in- 
tegration of the corresponding real-space overlap integrals. The relative sizes of these 
numerical contributions a posteriori justify the approximation of neglecting the color 
Coulomb and linear terms in 1=2 tttt scattering. 

The 1=2 TTTT scattering length has been calculated previously using many other theo- 
retical approaches. A summary of some of these predictions is given below. (We use a 
current value of = 93 MeV in Weinberg's PCAC formula a^=2 = -M^/IGtt/^I) 



J=2 



thy. 



' -0.053(7) [fm] LGT [7] 

-0.052 [fm] meson exchange p5| , pi 

-0.053 [fm] Roy Eqs. 

[ -0.063 [fm] PCAC |26 



(9) 



Although and other effective lagrangian parameters are normally taken from exper- 
iment, these parameters are of course calculable from quark-gluon forces. As an example, 
our result for Oq"^ yields the following expression for 



1 _ 2^ / 2^ 



m 



2 + 32 \^3l/2 



11^ +-(1 

(3^ 3 ' 



31/2 ^^ vr6 



(10) 



The dominant contribution comes from the 0{as/m?) OGE S ■ S term. 

Experimental determinations of the scattering length have yielded results which are 
larger than theoretical expectations; 

r -0.13(2) [fm] Losty et al. [21] 

"° expt. \ -0.24(2), -0.22 [fm] Prukop et a/, g 

We speculate that this discrepancy is due to the use of a simple 5 = kT,a + 0{kJ^) effective 
range formula in extrapolation. The difficulty of extrapolating experimental phase shifts 
to threshold has been stressed by Morgan and Pennington p7|,p8|. We advocate the use 
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of a "generalized specific lieat plot" of low energy phase shifts for this purpose . This 
plot takes into account the threshold behavior seen in Eq.(6), 

5^=^""^ = KE^ f{as,m,h,(3,x) (12) 

where / is a relatively slowly varying function of x = k^^'^/AjS'^. Thus the threshold 
behavior is approximately proportional to k-j^E^^ rather than just /Ctt, and since M^^ is quite 
small this leads to rapid variation near threshold and makes the linear-/c7r extrapolation 
inaccurate. We suggest removal of all this dependence by displaying 5o/(^7r-£'7r/^7r) versus 
/Ctt^- The intercept in this plot is the scattering length, and the slope at intercept implies 
the effective range. 

This generalized specific heat plot is shown in Fig. 5 for 1=2 vrvr scattering. An ex- 
trapolation of the moderate-A;^ data can now be seen to be much closer to the theoretical 
scattering lengths. The small-/c7r dependence of do/kj, was calculated by Donoghue 
in a chiral effective lagrangian, which gave the Weinberg result at fc^r = and an O(fc^) 
correction factor of (1 + /c^^/2m7r^). This is just the correction due to an overall factor of 
E'tt, so this predicts a zero slope in A;^^ for 5q/ {k^^ET^/MT^) at threshold. 

The Jiilich meson-exchange model which is dominated by t-channel p exchange 



in this process, also predicts rapid variation in 60/ kj^ near threshold. The prediction of 
this model for SQ/lkj^ET^/Mj^) |2^, shown in Fig.5, is rather similar to our quark model 
result. 



4- 1=2 vrvr Equivalent Potentials 

Low energy "phase shift equivalent" Gaussian 1=2 vrvr potentials, derived using the 
method of Mott and Massey as described in App.E, are given below. We quote 
separate Gaussians for the transfer and capture contributions from each of the three in- 
teractions, spin-spin contact, color Coulomb and linear hyperfine. However their predicted 
phase shift decays more slowly at large momentum, probably due to the use of power law 
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form factors in their vertices. 



"^ 327rl/2 m2 

I 29/2 ^ ^-|/32r2 
"^ 327ri/2 m2 

211/2 -6/32^2 

-31/253/2^1/2 e ^ 



2I/23I/2 

- ^1/2 e 8 



S/32r2 



29/23V2 b ^-f/32^2 
'73/2^1/2 /3 c 

2V235/2 g-A^2^2 



5* • S* (transfer) 
S ■ S (capture) 
Cou. (transfer) 
Cou. (capture) 
lin. (transfer) 
lin. (capture). 



(13) 



■53/2,11/2 /3 

In Ref. we derived 1=2 tttt potentials for the spin-spin contact interaction using the 
"locahty expansion" method of Ref. ||31]]. This gave an identical result for the spin-spin 
transfer potential, because this amplitude (before Bose symmetrization) is a function of 
t only. However for the capture diagrams the Mott-Massey approach used here gives a 
different potential, since it is constrained to reproduce the 0{k^) series expansion of the 
phase shift in Eq.(6), but the local approximation is not. The two capture potentials 
reproduce the scattering length, but the local approximation gives an incorrect effective 
range. 

The low energy Mott-Massey 1=2 tttt potential is shown in Fig. 6 for our standard 
quark model parameters as = 0.6, b = 0.18 GeV^ and m = 0.33 GeV. The spin-spin 
hyperfine contribution is dominant over the range shown. 



5. 1=2 TTTT Phase Shifts with L > 



The higher partial waves (L > 2) may be evaluated similarly. According to Eq.(^, 
these receive contributions only from the transfer diagrams. The Born-order D-wave phase 
shift with SHO wavefunctions is given by 



Xl=2 TTTT 
"2 



fc-E'TT ^2 I 



1 l-e-2a^ _|_ 1 l+e-2^ 



1 l-e-2^ 
3 a:3 



kE. 



71^2 

b 



32 X ' 3 X 

2 f^_i(-2x)-e-2- ^ 2 f^_j_(_2x)+e 
3^ 



+ 



32 



2 

3-5 



_3(-2x)-e- 
■ 7 



■ f, 1 (-2a;)-e-2^ 
\ 3 



1 ^2, 



1 {-2x)+e- 
_3 



^ 5-7 



L _3(-2x)-e- 



S-S 

Cou. 
lin. 

(14) 
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These three expressions are numerically rather small, and their phases are such that they 
approximately cancel; at M^^ = 1.5 GeV they are respectively —0.8°, +0.2° and +0.4°. 
To see this more clearly, the leading 0(/c7r^) behavior predicted by Eq.(14) is 

- ( - ^ ^ ^ ^ 4 ) t • ''^> 

and the three dimensionless combinations a^/^^/m^, and b/jS'^ are comparable in size. 
We have also evaluated this D-wave phase shift using Coulomb plus linear plus hyperfine 
wavefunctions. The result is shown in Fig. 4, and is numerically similar to the SHO D- 
wave, Eq.(14). 

In comparison the experimental D-waves reported by Durusoy et al. [|T9| and Hoogland 
et al. I^O] are ^ —3° at Mtttt = 1.5 GeV (see Fig. 4). (Losty et al |21[] report a rather larger 
but inconsistent low-energy D-wave.) This is clearly larger than our prediction, although 
the rather slow variation of the Durusoy et al. and Hoogland et al. D-waves with Mtttt 
may indicate a problem with the measurements; the expected threshold behavior of kj* 
is much more rapid than the observed energy dependence. Unfortunately the dispersion 
relations represented by the Roy equations have technical difficulties with determining D- 
and higher waves [^]. They do however lead to predictions of a positive D-wave close 
to threshold, which is not evident in the data. The D-wave may well have important 
meson exchange contributions, since this type of model can accommodate the reported 



experimental phase shift p2|. 



C. PsV Scattering 

1. 1=2 Tip T-matrix 

For simplicitly we will initially quote results only for 1=2 vrp. The other isospin 
channels are simply related by flavor factors, which we will discuss subsequently. We 
assume identical spatial wavefunctions, so only the p spin degree of freedom and difference 
in phase space distinguish this case from tttt. Summing the individual contributions in 

14 



App.B with the appropriate flavor and color factors and the (—1) signature phase, and 
using the PsV spin matrix elements of App.C part C3, we find for the 1—2 irp T-matrix 



T 



1=2 irp 



fi 



39/2 



-A 2/3/32 



-A 2/2/32 



Sp-i{AxC) 



Sp-i{AxC) 



+ 



m^/S^ I 32-5 2-2 



-A 2/2/32 



Sp • Q— Sp ■ Q— ttQ- 



+^( + Jr^ h 1(^76/9^) 1 e-^'/'f^' 



m- 



!/J2 I 39/2 . 5 2'2 



— * — * — * — * 2 — *2 



+ 



Sp-CSp-C--C' 



(16) 

The individual contributions in this result are respectively OGE spin-spin; OGE color 
Coulomb; linear confinement; OGE spin-orbit; linear scalar confinement spin-orbit; OGE 
tensor (transfer diagrams); and OGE tensor (capture diagrams). In all these we list 

— * — * — * 

transfer followed by capture contributions. Sp is the p meson spin vector, A and C 
are the initial and final tt momenta, B — —A and D — —C are the initial and final p 

— * — * — * 

momenta, and Q± = C ± A as in tttt. Since this result was derived in the CM frame, 
1^41 = \B\ = \C\ = \D\. This T/j evidently describes irp spin-orbit and tensor interactions. 
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in addition to spin-independent scattering. It is interesting that there is a one-to-one 
mapping between the quark-quark spin-orbit and tensor interactions and these vrp spin- 
orbit and tensor terms. This simple result need not be true in general; a given spin- 
dependent interaction at the quark level may give rise to a different type of hadron- 
hadron interaction. As an example, a mapping of a tensor nucleon-nucleon force into a 



nucleon-nucleus spin-orbit interaction was discussed by Stancu, Brink and Flocard. ||33 

To evaluate phase shifts and inelasticities it is convenient to calculate the matrix 
element of our vrp T- matrix Eq. (|l6|) between general \jls) states, which gives the reduced 
matrix element 

Tin = {3m,l's\T\jm,ls) = 
J2{jm\l'ix',ls:){jm\lix,ls.) jj dil' dil Y;^,in') (l4|T;,(fi', ll^(fi) , (17) 

mm' 

as discussed in App.D. This is a straightforward exercise, although integrals of special 
functions and a careful angular analysis of the spin-orbit and tensor terms are required; 
the details are discussed in Apps.G and H. This matrix element is diagonal in / except for 
the tensor interaction, which has both diagonal and off-diagonal (transfer) and fully off- 
diagonal (capture) contributions. The /-diagonal results, again showing transfer diagram 
contributions followed by capture, are 



nicou. = ^ (-Sl,e.en ^ -^E (^) + |^ /i,|(2x/3)') C''^ (19) 



^4 y(>l,even g 



mi^. = ^ (Sl,e.en - :F%{x) - 6,0 /-i,|(2x/3) 6"^^ (20) 
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OGE L-S 



TV? 

25 



(L-S) 
1 



28 



32 (2/+1) 



X 



i^m^ - (x)) -S,,^x /3,|(2./3) e-- (21) 



lin. L-S 



71% 
l4 



(L-S) 



2x 



rpj I transfer 
^ll\OGE T 



(T) 



-1) 



(22) 



X 



-^tl\^) + 77^4'^^) + J^^!.7J:^ o^ ^f/H^)l e-^^ (23) 



33-5 V(2^ + l) ^ (2Z + 3) f'i^ ' (2Z + l)(2Z + 3) U 

and the off-diagonal tensor matrix elements are 



rpj I transfer ^ C^s 



(^ZJ-l^/'J+l + '^/J+l'^rj-l) 

xf4^r^(x) + 2^^;(x)+nT'(a;)) e-^^ (24) 



I capture ^ '-'^s 



32-5 (2j + l) V 5.5 



(^ji [Si2Si'o + Siq5i'2) g^^ya . 5 ^ /|,|(22^/3) e 



-2i; 



(25) 



In these formulas ii{x) is a modified spherical Bessel function, the tensor (T) matrix 
element between \j,l,s = 1) np states is 

- 1 J =1+1 

(T) = J -(21 + 3)/l j = I (26) 

.(/+l)(2/ + 3)//(2/-l) j = l-l , 



and the integral 



J ^ dn iFi(a; c; a;(l + /x)) 



(27) 



is evaluated in App.G. 
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2. 1=2 7rp S-wave Phase Shifts 



In S-wave to S-wave scattering the spin-orbit and tensor np T-matrix contributions 
vanish, and we are left with color Coulomb, linear and spin-spin contributions, just as 
in the 1—2 tttt case. The 1=2 7rp S-wave phase shifts that result from these interactions, 
again using Eq.(D17), are 



.7=2 





as 




m? 


kE^Ep 


ds 


VS 




kE„Ep 


b 


VS 





(-SKi-'-")- 



26 _-4a:/3 
WPi c 



s-s 



3, - (f,, . (-2x) - e-^-) - ^ f,, 3 (-2X/3) e-^-/3) Cou. (28) 



32 



2f. 1 



-2x) 



) 



22 r 



3172 ^2 i 



-2a;/3) e 



~42;/3 



lin. 



T*2 



where -^/i = {Et^ + Ep), and again x = A /4/3^. In Fig. 6 we show these individual 
components and the total S-wave phase shift with our standard quark model parameter 
set and meson masses (used throughout) of M^r = 0.138 GeV and Mp — 0.77 GeV. The 
forces considered here evidently lead to strong repulsion in the 1=2 irp channel. 



3. 1=2 irp Phase Shifts with L > 

The spin-orbit and tensor terms in Eqs. (21-23) all contribute to Z > 7rp scattering, 
and there is also an odd-/, j-independent term due to the OGE spin-spin interaction 
in Eq.(18), which is not symmetric under p —p. The color Coulomb and linear 
confinement spin-independent terms, Eqs. (19,20), contribute only to even I. 

Adding the various diagonal matrix elements of Eqs. (18-23) and using Eq.(D17) gives 
phase shifts for each ^Lj partial wave. In Fig. 8 we show results for all P-wave channels 
and for J=L±1 in D- and F-wave. Note that there is a large, inverted spin-orbit force in 
the P-wave, so the ^Pq phase shift is widely separated from ^P2, and has an even larger 
maximum phase shift than the S-wave. The higher-L channels show decreasing phase 
shifts with increasing L, as expected for short-ranged quark-gluon forces. 

The relative importance of the individual contributions to the spin-dependent force 
is of considerable interest. In Fig. 9 we show the various spin-dependent contributions 
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to the 1=2 vrp phase shift. The largest contribution arises from OGE spin-orbit, 
in particular from the transfer diagrams. The OGE and confinement spin-orbit capture 
diagrams give smaller contributions of the same sign. Finally, the confinement spin- 
orbit transfer diagrams have a sign opposite to all these and reduces the total spin-orbit 
force somewhat. This dominance of the PsV spin-orbit by OGE is an interesting result. 



especially since Mukhopadhyay and Pirner |3J| found the opposite result in KN. In that 
system they concluded that confinement, not OGE, makes the largest contribution to the 
spin-orbit force. The OGE tensor in 1=2 '^P2 vrp is weakly repulsive; it makes a much 
larger contribution to '^Pi and '^Pq, where the tensor matrix element is respectively —5 
and 10 times as large. The OGE tensor is evident in Fig.8, in the departure of the ratio 
(3p^ — •^Pi):('^Pi — ^Pg) from the pure spin-orbit value of 2:1 at higher energies. 

There is also an off-diagonal coupling due to the OGE tensor terms, given by 
Eqs. (24,25), but we have neglected this in calculating phase shifts because we find that it is 
numerically a small effect. The largest coupling at low energies is ^Si ^ '^Di, which leads 
to an inelasticity of only rjsD = 0.97 by M^^p = 3.0 GeV (calculated using Eqs.(D18-D20)). 



4- 1=2 np P-wave Spin-Orbit Potentials 

We may determine low energy Gaussian equivalent vrp potentials from the phase shifts, 
as discussed in App.E. The most interesting potential phenomenologically is the spin-orbit 
one, since the origin of the spin-orbit interaction in the NN system is a long-standing and 
still poorly understood problem. In particular we derived Gaussian potentials correspond- 
ing to the P-wave phase shifts due to the OGE and linear scalar confinement spin-orbit 
interactions, using Eq.(E4) of App.E. The results for the transfer and capture contribu- 
tions to these potentials are 
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275/2^1/2 rn? 
55/2 aj^ 



- 3275/2^1/2 ^;;2- (L • S) e 7 /3 OGE (transfer) 



■3972^ — (L-S) e t^''-' 



i^2^2 



(29) 



55/2 

'21/275/2^1/2 ^2 



^ (L ■ S) e- 



OGE (capture) 
lin. (transfer) 
lin. (capture) . 

The OGE, linear and total spin-orbit potentials for the '^P2 wave of 1=2 vrp are shown 
in Fig. 10 for our standard parameter set. The largest contribution to the vrp spin-orbit 
force comes from OGE transfer diagrams; the linear confinement spin-orbit from the 
transfer diagrams is about half as large and opposite in sign, and the two capture diagram 
contributions are much smaller. Since the confinement capture and transfer diagrams have 
opposite signs, the net result is dominance of the PsV spin-orbit by the OGE contribution. 



D. Scattering in J^^" -Exotic PsV Channels. 

TABLE I. J^"^" Exotic States in PsV. 



Channel Exotic Quantum Numbers 



Meson Pair 


hot 


s 


p 


D 


F 


G 


Tip 


0, 2 






2+- 




4+- 




1 




1-+ 




3-+ 




ITU), rjp 


1 




0— 


2+- 




4+- 


rjuj 









2+- 




4+- 



The recent evidence for J'^'^"-exotic resonances 7ri(1400) and 7ri(1600) I^S] has made 
the study of scattering amplitudes in exotic channels especially interesting. The surpris- 
ingly low mass of the 7ri(1400) in particular has led to suggestions that it might not be 



a "hybrid" gluonic excitation, since these are expected at ~ 1.8 — 2.0 GeV PBI. Another 
possibility that the 7ri(1400) is a "multiquark" , perhaps a meson-meson bound state in 
a very attractive channel. We can test the plausibility of this type of assignment by 
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calculating meson-meson scattering amplitudes in the various exotic channels. 

The exotic channels accessible to the hghtest nonstrange PsV meson pairs are listed 
in Table I. (We do not tabulate light PsPs exotic amplitudes because they are zero in this 
model. The PsPs exotic channels are odd-Z ttt], ttt]' and 7777', whereas the quark interchange 
model PsPs scattering amplitudes are even-/, assuming identical spatial wavef unctions.) 
We generally expect the largest scattering amplitudes to be in the lower partial waves. 
In PsV the P-wave has the first exotics, which are J^*-^" = (all channels except 1=1 
np) and J^*^" — 1"^ (1=1 np only). Calculation of these scattering amphtudes simply 
requires changing the external qq flavor states attached to the Feynman diagrams of Fig.l. 
The results relative to the 1—2 np case treated in the paper are summarized in Table II. 
TABLE II. Overall flavor factors in diagonal PsV scattering. 



Channel Relative Amplitude 



irp 


1=2 


+1 




1 










-1/2 


TTLO 


1 


+1/2 


VP 


1 


+1/4 







+1/4 



Inspection of the tables shows that the largest exotic scattering amphtude should be 
in the 1=2 vrp P-wave. The elastic phase shift in this channel is the ^Pq curve in 
Fig. 8. The large negative phase shift shows that this is a strongly repulsive channel; the 
maximum phase shift is predicted to be a quite large fa —50° at M^p ps; 3.1 GeV, which 
exceeds even the S-wave phase shift maximum. The largest attractive exotic phase shift 
we have found in PsV is the 1=0 partner, which is —1/2 of 1=2, giving a maximum phase 
shift of ~ +25° at the same mass. We do not find sufficient attraction to form a meson- 
meson "molecular" bound state in any of these nonstrange J^'-^"-exotic PsV channels. 
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The 7] channels are relatively weak because only the nn part of the t] contributes to these 
diagonal scattering amplitudes; the ss component leads to open-strange final states (K*K 
for example) after quark line interchange. 

Regarding candidate exotic resonances, there have been speculations that the determi- 
nation of the mass and width of the exotic candidate 7ri(1400) may have been compromised 
by inelastic rescattering effects |^ , analogous to the "Deck effect" proposed as a nonres- 



onant explanation of the ai(1260). For example, crossing the nbi threshold at ~ 1.4 GeV 
in the process vrp ^ vrfoi tti] might mimic resonant phase motion if this process has 
a rapidly varying inelasticity. We can test this and other nonresonant possibilities by 
calculating the elementary 2-^2 scattering amplitudes using our quark model approach. 

Some important results follow from simple flavor factors. Note in particular that the 
nonresonant scattering amplitude np ^ np vanishes in any 1=1 channel, including the tti 
exotic one. This is a general result whenever the quark line diagram of Fig.l dominates; 
clearly a pair of oppositely charged, nonstrange qq mesons A^B^ cannot scatter into 
another charged pair C^D~ under quark interchange. A comparison with isospin matrix 
elements shows that this implies that scattering of any two qq isovectors in 1=1 vanishes. 
This isospin selection rule eliminates two subprocesses discussed by Donnachie and Page 



| 37| as Deck effect backgrounds that might shift a higher- mass exotic resonance to an 
apparent 7ri(1400), vrp irbi —>■ ttt] and 7rp ^ vrp ^ Tcrj. 

Independent of any scattering model, one should note that the coupling np nbi is 
probably small because of the strong VES experimental limit (reported by V.Dorofeev 
@)of 

5(7r2(1670) ^ 7r6i) < 0.19% (2a c.l.) . (30) 

Since tt2 vrp is a large mode {B = 31(4)% [^), if vrp nbi rescattering were important 
we would also expect to observe a large 112 vrfoi branching fraction. 

We also expect the flnal background process suggested by Donnachie and Page (vrp — > 
711] 711]) to be small, because the direct time ordering 7n] 711] vanishes in P-wave in 
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this model. Finally, the rescattering process they propose, 7ri(1600) nbi rjn, does 
not vanish in the quark interchange model, although the required AZ = 1 and suppressed 
1] flavor factor may nonetheless make this a relatively weak amplitude. A calculation of 
this and related scattering amplitudes is planned for a future publication. 

E. Experimental prospects for measuring PsV phase shifts. 

Although there is little experimental information about PsV interactions at present, 
these phase shifts actually are experimentally accessible in existing data, for example as 
relative FSI phases in the D and S amplitudes in bi —>■ uj-k. These are usually, and in- 
correctly, taken to be relatively real amplitudes. The relative phase including the FSI is 
D/S = |D/S| ■ S^o-^s) [|8|^ and is observable for example reduction in the strength 
of the SD cross term in the tiuj angular distribution by cos{5s — ^u)- Since this method 
requires individual measurements of the S^, and SD cross term in the angular distri- 
bution, it should be applicable to cases such as 61 (1230) ttu and 61 (1600) —>■ ttuj where 
S and D are of comparable magnitude |^. The 630 = — ^d phases we predict at these 
masses (which are calculated from +1/2 times the ^Si and '^Di 1=2 np phases in Figs. 7,8) 
are S^'^iM^^ = 1.23 GeV) = -14.° and S^'^iM^^ = 1.60 GeV) = -17.°. 

This proposed technique is similar to that used in Ke4 decays |^, in which the low 
energy 1=0 tttt S-wave phase shift is actually observed as the difference between the 1=0 
S-wave and 1=1 P-wave tttt FSI phases. 

III. SUMMARY AND CONCLUSIONS 

In this paper we have derived meson-meson scattering amplitudes, including spin- 
dependent forces, from a calculation of the Born-order matrix element of the quark- 
quark interaction between two-meson states. Since qq annihilation is not included in 
the model, it describes scattering that does not involve coupling to s-channel resonances. 
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This includes for example 1—2 and the nonresonant backgrounds in all channels, including 
exotic J^'^". 

We considered the cases of PsPs and PsV scattering, and derived the scattering ampli- 
tudes in all j,l,s channels for these cases. The parameters of the model were previously 
fixed by quark model studies of hadron spectroscopy. Where possible we have compared 
the results to experiment. 

In 1—2 TTTT (the best studied PsPs case) the results were shown to be in reasonable 
agreement with experiment in S-wave scattering, and an extremum predicted near M^^r = 
1.5 GeV is supported by the data. Rapid variation of ^o/^tt is predicted near threshold, 
which may reconcile theoretical expectations of a small scattering length with larger 
reported experimental values based on extrapolation in kj^. The experimental D-wave, 
although quite small, is clearly larger than the model predicts. 

The PsV system is a convenient theoretical laboratory for studying spin-dependent 
forces, since it can accommodate both spin-orbit and tensor interactions, and is simpler 
than KN or NN. We derived analytical results for these spin-dependent PsV interactions 
(T-matrices and phase shifts) given SHO wavefunctions and the standard spin-dependent 
quark model forces. The quark-quark spin-orbit and tensor forces map directly into spin- 
orbit and tensor PsV interactions. We find that the OGE spin-orbit force in the PsV 
system is quite large in P-wave, and so is expected to be large in many other hadron- 
hadron systems as well. 

There is no PsV phase shift data at present. We noted however that PsV phase shifts 
actually can be measured in multiamplitude resonance decays to PsV final states, so it 
should be possible to test theoretical predictions for PsV scattering amplitudes in future 
experimental studies. 

Our predictions for scattering in J^*""-exotic channels are of current interest because 
the reported exotics might be complicated by large and rapidly varying nonresonant in- 
elasticities. One speculation is that the 7ri(1400) parameters might be strongly affected 
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by the opening of inelastic couplings to the irbi channel. In our model (and in any qq 
constituent interchange model) several of these nonresonant processes can be rejected as 
significant complications because of vanishing flavor factors. 

In future we plan to extend our calculations to other exotic meson-meson channels, 
such as S+P, to test whether strong attractive interactions are predicted that might 
support "multiquark exotics" such as S+P molecules. We also plan to apply the current 
approach to the study of spin-dependent interactions in other hadronic systems, including 
KN, NN and light hadron + charmonium systems. 
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APPENDIX A: QUARK-LEVEL T-MATRICES AND WAVEFUNCTIONS 

The various contributions to the quark-quark Tf^ (with color factors of T^T" removed) 

are 



Tfi{q,Pi,P2, 



8nas 
4-7ras 



I 67rb 

I Aiwas 
3inb 



Si ■ S2 



Si-{qx{^^-p2))+S2-{qx{pi- 

— * — # 

Si - {qx pi) - S2- {qx P2) 



Si-qS2-q-kq^Si-S2 



OGE spin-spin 
OGE color Coulomb 
linear conft. 
OGE spin-orbit 
linear spin-orbit 
OGE tensor. 

(Al) 



The standard qq quark model Gaussian wavefunction is given by 

1 



7^3/4^3/2 



where in general 



Prel 



rUqPq - rUqPq 



{iriq + mq)/2 

and for our special case of equal quark and antiquark masses 



(A2) 



(A3) 



Prel =Pq-Pq 



(A4) 
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APPENDIX B: EXPLICIT OVERLAP INTEGRALS 



Bl. Results Included 



In this appendix we give the exphcit meson-meson T-matrix elements that follow from 
the overlap integrals Eqs.([I|-^) with Gaussian wavefunctions and the various quark T- 
matrix elements. The OGE spin-spin hyperfine, color Coulomb and linear confinement 
results were derived previously [|l^ . For completeness we quote the formulas here, as well 
as giving the new spin-orbit and tensor results. The multiplicative diagram-dependent 
color and flavor factors and the signature phase (which is (—1) for these meson-meson 
scattering diagrams) are not included in the results given below. These formulas abbre- 
viate the confluent hypergeometric function iFi(a; c; x) as fa,c(a^), and Q± = {C ± A). 



B2. OGE Spin-Spin Hyperfine Contribution 



These simple contact matrix elements were evaluated previously, for example in Ref. |^ 
(in an equivalent form, but incorporating color factors and the signature phase, as Eqs. (TI- 
TS) of that reference). The results are 



2^ vra. 



3 m2 



Si ■ 5*2 



2^ vra 

35/2 ^2 

rr^.{C2) _ rp (CI) 



£ g-A 2/3/32 



Si ■ S2 



fi 



(Bl) 
(B2) 

(B3) 
(B4) 



B3. OGE Golor Goulomb Gontribution 

The contribution of the OGE color Coulomb interaction to the meson-meson T-matrix 
follows from the evaluation of the integrals Eqs. (0-|5) with the second quark-quark Tfi in 
Eq. ([Mp . The results are 
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T/^'' = + 2^ ^ fi,|(Q!/8/3^) e-^'/'^' (B5) 

T;.(^2)=T;,(T^)((7--C7) (B6) 

T//^^^ = + |i ^ fit (^^V6/?^) e-^'/'^' (B7) 

T,/^^) = T,.(ci) . (B8) 



B4- Linear Confinement Contribution 

The linear confinement integrals were carried out in coordinate space, since the Fourier 
transform of the linear potential is singular. The results are 

T/.(T^) = -6 ^ L,,iQl/8P') e-^'/'^' (B9) 

P 2 2 

Tf,(^^)=Tf,m^C^-C) (BIO) 
T^,(ci) = _ 33/2 ^ f , 3 (12/6/32) e-^^'/2/5^ (Bll) 



One may also obtain these results using the momentum space integrals Eqs.(|l]-^, but 
the quark-quark Tji in Eq. (|AlD must include a long-distance regularization in the 
intermediate stages of the integration. The final result is well defined due to the Gaus- 
sian damping provided by the wavef unctions. We have checked both the linear and color 
Coulomb Tfi results by comparing the expressions Eqs.( p5| - P12|) with Monte Carlo eval- 
uations of the corresponding real-space overlap integrals. 

B5. OGE Spin-Orbit Contribution 

The four OGE spin-orbit overlap integrals can be evaluated similarly using the fourth 
quark-quark Tfi in Eq.([Al|), which gives 
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T..(T1) 



m?[3'^ 2'2 



f3 5(g!/8/32) e-^''^^' 



{Si + ^2) ■ i{A X C) 



35/2 ^2^2 5-2 
J- ft — J- fi 



(Si - S2) ■ t{A X C) 



(B13) 
(B14) 

(B15) 
(B16) 



B6. Scalar Confinement Spin- Orbit Contribution 



The matrix elements Eqs.(|l|-^), of the scalar confinement spin-orbit interaction in 



Eq.i 



are 



{Si + S2) ■ i{A X C) 



nb 



2 ■ 31/2 m^pi 2.2 

T,.(C2) _ _ rr. (Cl) 



{Si-S2)-i{AxC) 



(BIT) 
(B18) 

(B19) 
(B20) 



B7. OGE Tensor Contribution 



Finally, for the OGE tensor Tfi (the last entry in Eq. ([Al|) we find 

1 



5 m?j3^ 2'2 ' 



Si-Q.S2-Q---QtSi-S2 



r//^2) = r;,(T^)((7--(7) 



2^ 



' 35/2.5 ^2^2 

T-^.(C2) _ T-^.(Cl) 
-L fi — -^fl 



Si- AS2- A- -A^ Si- S2 
3 



(B21) 
(B22) 
(B23) 
(B24) 



The tensor matrix elements in the capture diagrams, Eqs.( p23| , p2l) , are the only 
cases in which we have found a post-prior discrepancy in these PsPs and PsV scattering 
amplitudes; the post forms of these matrix elements involve a tensor in C rather than A, 
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T^. (CI, post) _ I 



35/2 . 5 ^2^2 2-2 



f|r(A76/3^)e 



2n ^-12/2/32 



o\ ' 02 ' ^ ~ 01*02 
o 



(B25) 



These capture tensor terms vanish in the PsPs channel. They do make a small, off-diagonal 
contribution to PsV scattering, alheit only in the ^Si ^Di amplitude. 
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APPENDIX C: MAPPING QUARK SPINS INTO HADRON SPINS. 



CI. Spin Matrix Elements 

In these scattering amplitude calculations the matrix elements of spin-dependent quark 
interactions (the spin-spin, spin-orbit and tensor forces) involve matrix elements of linear 
and bilinear quark spin operators. Since the quark spins are not directly observed, it 
is useful to replace them by the spins of the external hadrons. This appendix gives the 
(diagram dependent) mapping from quark spins to hadron spins in the PsPs and PsV 
cases considered in this paper. 

The spin matrix elements we require are /, -S'(l)*, 5'(2)' and S{lYS{2y between general 
initial and final PsPs and PsV spin states. Our convention for the diagrams (Fig.l) is 
that mesons A and C are always Ps {e.g. tt), and B and D are Ps or V {e.g. tt or p). 



C2. PsPs 

First, in PsPs scattering there are no external meson spins, so the quark spin matrix 
elements are proportional to geometrical tensors such as The matrix elements by 
diagram are 

1 

{PsPs\ I \PsPs) = +- all diagrams (CI) 



{PsPs\ S{iy \PsPs) = {PsPs\ S{2y \PsPs) = all diagrams (C2) 



~I~Q ^ Tl, T2, Tlgj^j^j^, T2,gyffiffi 
{PsPs\S{iyS{2y\PsPs) = { ^ hence (C3) 

8 ^ C2, dsymrm ^2gyjjijji 

I 3 rpi -yi -yo 

{PsPs\S{l)-S{2yPsPs) = { I . (C4) 

g CI, C2, C\symmt ^2gyjjijji 

Note that the spin-orbit and tensor terms vanish identically in PsPs scattering; this follows 



from applying Eqs.(C2,C3) to Eq.(Al). 
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C3. PsV 



In PsV scattering the vector {e.g. p) meson spin Sp provides an additional degree of 
freedom, and the hnear and quadratic quark spin matrix elements can be expressed in 
terms of the p spin matrix elements {pf\S'^p\pi) and The mapping of quark 

to meson spins is 

((7rp)/|7|(7rp)i) = +^ (p/|7|pi) all diagrams (C5) 

(C6) 

~l~4 \Pf\^p\Pi) T2, TUgyjjijji, C2, Cisymm 

{(7rp)f\S(2y\(7rp),) = +^ {pf\S;\p,) all diagrams (C7) 



{{np)f\S{iyS{2y\{np),) = 



f +I6^^pf\l\p,) + l^e^^>^ (p/|5>.) - i {Pf\S;Sj,\p,) Tl, T2,,„^ 

+ 1 + I z e^^'^ (p/I^JIa) T2, T1,,^^ 

-I a) - I i e'^^ {pf\Sl\p,) CI, C2,,^^ 

[ -|<5^^(P;|/|a) - i^e^^-^ + \ {Pf\SlSl\p,) C2, Cl,,„^ 



(C8) 
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APPENDIX D: PHASE SHIFTS AND INELASTICITIES FROM THE 

T-MATRIX 

Since total angular momentum is conserved, the T-matrix is block diagonal in a total 
angular momentum basis, and can be written as 



T = Y^ \jm) {jm\ . (Dl) 

jm 

The coefficients {Tj} can be determined by evaluating the matrix element 

T, = {jm\ T \jm) . (D2) 

In the special case of spinless scattering these basis states are eigenstates of orbital angular 
momentum 

\im) = Jdn Y^n) \n) , (D3) 

so the T-matrix is given by 

T = ^ |/m) Ti {lm\ =J2jJdQdQ' \Q') Yi,^{Q') Ti Y;^{Q) {Q\ . (D4) 

Im Im 

The T-matrix can also be written in terms of the scattering amplitude Tfi{Q, Q') between 
momentum eigenstates, 

T = JJ dn dVL' Tfi{n,n') {n\ , (ds) 

so Tfi{Q, Q') and Ti are related by 

I m / 

and hence 

Ti= {lm\T\lm) = jj dQdQ' Y;^{n') Tf,{n,n') Yim{n) (D7) 

where /i = cos 6nn' . A more familiar quantum- mechanical result follows from fixing the 
incident direction fl = z in Eq. (|D6| ) and integrating over final angles with a -Pz(/i) weight, 
which gives 
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Ti = 2ti j'dfiPi{fi)Tfi{z,n') . 
Since we define Tfi{Q,fl') by 

{Cn', D_n' I S \An, B_n ) = 5j,-i {2Txf 5^*\A + B - C - D) Tf, 



(D8) 



(D9) 



it is related to the Lorentz invariant 2-^2 scattering amplitude Ai defined by the PDG 
(their Eq.(35.8)) by 

M 



and hence to the CM scattering amplitude f{k,9) (their Eq.(35.48)) by 



(DIO) 



fi 



(Dll) 



(Here k is the CM momentum of any particle, k 
wave expansion of f{k,9), Eq.(35.44), 



\A\ 



\B\ = \C\ = \D\.) The partial 



fik,e) = J2i2l + l)aiPiifi) 



(D12) 



and the relation between a diagonal partial wave amplitude a; and the phase shift 



ai 



2i 



(D13) 



allow us to determine 6i from Tji . For purely elastic scattering, and assuming small phase 
shifts so that ai ^ 6i, Eqs.( P6| - P13| ) give 

1 kEAEs 



in^ ^/s 



Ti . 



(D14) 



This relation was used previously to determine for example Kvr [11] and Z-diagonal KN 



1 12 1 elastic phase shifts. In the case of elastic scattering of identical bosons, such as 1=2 



TTvr, there is an additional factor of two for identical particles |38], so the relation between 
the Born-order phase shift and the T-matrix element becomes 

kEA 



5i 



ident. 



(D15) 
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Since the angular integral in Eg. ( p 141) is proportional to the amplitude {lm\T \ lm) , 
we may also write the Born-order elastic phase shift (for distinguishable particles) directly 
in terms of the T-matrix, 

1 kEAEs 



87r2 



{lm\ T \lm) 



(D16) 



This formula has a straightforward generalization to the case of external hadrons with 
spin, which we use to evaluate irp phase shifts and inelasticities. In the case of an /- and 
s-diagonal interaction this is 



Stt^ 



(D17) 



This is adequate for diagonal forces such as our vrp spin-orbit interactions. The ten- 
sor force however is not /-diagonal, so for this case we must introduce a more general 
parametrization. Since the tensor interaction couples /, /' channel pairs which have the 
same j but differ by |/ — /'| = 2, it leads to 2 x 2 iS-matrices. These can be parametrized 



as 



(D18) 



In our calculation, both of the phase shifts and the inelasticity ew = y 1 — rifi, are 
0{Hj), so to this order we can relate these linearly to the matrix elements of Hj. The 
Born order phase shift formula (D17) remains valid for both channels, and the inelasticity 
is 



^11' 



1 kEAEs 



47r2 ^ 



{jm; l's \ T \jm; Is) . 



(D19) 



The overall phase of eur is dependent on the state normalizations, but the familiar 



Vw 



^ - p2 



(D20) 



IS unique. 



35 



APPENDIX E: POTENTIALS FROM THE T-MATRIX 



Potentials provide a very useful representation of hadron-hadron interactions. They 
have a clear intuitive meaning, and can easily be used in the nonrelativistic Schrodinger 
equation in searches for bound states or in coupled channel problems. 

Unfortunately, one may define hadron-hadron potentials in many different ways. Ide- 
ally they should reproduce phase shifts or T-matrix elements, at least in the low energy 
limit. The assumption of a unique, purely local potential is in general overly restrictive, 
as it leads to a scattering amplitude that is a function of t only. In general we find 2 — 2 
scattering amplitudes that depend on both s and t. One approach to this problem is to 
introduce nonlocal "gradient" corrections to the potential |3l[] , which can be expressed for 
example as V{r) L ■ S terms; this approach leads to the familiar Breit-Fermi Hamiltonian 
for one-photon and one-gluon exchange, and was used to define hadron-hadron potentials 
in our previous work |^,|13 . 



Alternatively one may project the scattering amplitude onto a given angular channel 
/ so that only s dependence remains, and find a local /-wave potential that describes the 
scattering in that channel. This definition of potentials was discussed by Mott and Massey 
poll , and is equivalent to the definition we shall use here. This approach was previously 
used by Swanson |]Tn[ to define meson-meson potentials from scattering amplitudes. 

The quantum mechanical relation between the phase shift 6i{k) and the radial wave- 
function Ri{r) in potential scattering of a mass fi particle (which becomes the reduced 
mass below) in first Born approximation is 

6i = -2fik / r^dr V{r)ji{kr)^ . (El) 
Jo 

Since our T-matrix elements implicitly determine the elastic scattering phase shifts, for 
example the / = 2 tttt S-wave in Eq.(6), we can invert this formula for each / to determine 
the corresponding Z-wave local potential V/(r). In practice we find that our phase shifts 
are sufficiently "hard" at high energies to require singular potentials; this is presumably 
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an artifact of our approximations, such as assuming a contact spin-spin interaction. For 
this reason we do not completely invert the phase shift relation (El), and instead simply 
fit an assumed Gaussian form 

V{r) = Vg e-'^l'l (E2) 

to our theoretical low energy phase shift. The two Gaussian parameters are determined 
from the and 0{k^^ terms in the expansion of the phase shift near threshold, which 
in the S-wave case are equated to 

Ihn b,{k) = i^Vgrl k(l-lrie + 0{k')) . (E3) 

The generalization to higher Z, also using Eqs.(El), is 

Ita 5,{k) = -^-^fiV,,-! 1^ (l - 1 r^^e + 0(k*)) . (E4) 

In determining the Gaussian potentials that correspond to our derived phase shifts 
such as Eqs.(6,14,28), we set the external factors of and Ep equal to m^r and nip before 
expanding in k. This corresponds to using nonrelativistic phase space and nonrelativistic 
external hadron line normalizations in our T-matrix calculations, which we assume is the 
appropriate choice for the derivation of a nonrelativistic equivalent potential. 
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APPENDIX F: THE POST-PRIOR DISCREPANCY 



The "post-prior discrepancy" is a familiar problem in rearrangement collisions; the 
diagrams of Fig.l treat the scattering as due to an interaction between the initial hadrons 
A and B ("prior"), but we could equally well have written the interaction between the 
two final hadrons C and D ("post"). The post T-matrices may be obtained from the prior 
ones by exchanging the initial and final mesons and transforming the momenta A ^ C, 
C A and q —q. Thus for example the post CI T-matrix is 

j3„ j3 



Tf,{q,p,-p + A) <l>A{2p-q-A)<^B{2p + q-A-2C) . (Fl) 



y(CLprior)^^^ ^ CD) 

args A^C ,C^A,q^ — q 



(F2) 



One may show that the post and prior results for the scattering amplitude are equal 
provided that the external wavefunctions are eigenf unctions of the Hamiltonian ||41|| . 
Swanson ||10[ shows an example of convergence of post and prior results for meson-meson 
scattering amplitudes derived from quark Born diagrams as the external wavefunctions 
approach exact Hamiltonian eigenstates. Of course the Gaussian wavefunctions we use 
to derive our analytical results are only approximations to the eigenf unctions of the full 
OGE plus linear Hamiltonian, so in general we can expect to find a post-prior discrepancy. 
In this study we actually find such a discrepancy only in part of the tensor interaction in 
PsV scattering, for which we take the mean of the two results, 

Tf,{AB ^CD) = ^ (Tf°'\AB CD) + Tjr'{AB ^ CD)'j . (F3) 
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APPENDIX G: Y lm EXPANSIONS AND RELATED INTEGRALS 



It is useful in the partial wave decomposition of scattering amplitudes to expand 
functions of the sum and difference momentum transfers Q± = C ± A (here \A\ = \C\) in 
spherical harmonics, 

fiQl) = E fM') E Y;^{Qc)YeU^A) , (Gl) 

i rri 

fiQ-) = E(-l)'/^(^') E YLi^c)YeminA) . (G2) 

e m 

This expansion may be inverted to determine the coefficient functions {/^(A^)}, 

f,{A^) = 27T jj^i P,(/i) f{Ql) 

= 27r jjiiP,{ii) f{2A\l + ii)) =2ni-lY IjfiP.ifi) f{2A\l-fi)) . (G3) 

Many of the scattering amplitudes derived in this paper are proportional to confluent 
hypergeometric functions in Q'^ or , and their partial wave decomposition requires the 
integral of a Legendre polynomial times a shifted confluent hypergeometric function. This 
integral is given by (abbreviating iFi(a; c; x) as fa,c{x)) 

' U) (a)-„-l {fa-m-l,c-^-l{2x) + (-1)^+-+!) 



•^iS(^) = /V Pi{^^) /a,c(x(l + /i)) = E <^ 
"' — 1 n 



m / ^ m+l 

m=0 



(G4) 

where 

2™m! (£-m)! ^ ^ 



and the Pochhammer symbol of negative index is 

, , Via — n) 1 



(G6) 



r(a) nLi(«-A:)' 
In numerical evaluations it is often useful to transform the confluent hypergeometric 
function in Eq.( |G4D to a more rapidly converging negative- argument form, using fa^d'^x) = 
fc—a,c{ 2x) e 
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APPENDIX H: SPIN-ORBIT AND TENSOR MATRIX ELEMENTS 



As noted in App.D, when evaluating phase shifts and inelasticities it is useful to 
determine T-matrix elements between \ jls) states. In the PsV system there are spin-orbit 
and tensor contributions to the T-matrix, and determination of the j, I, s-basis matrix 
elements of these terms is a complicated problem in angular analysis. Here we show how 
these matrix elements may be evaluated. 

First consider the spin-orbit terms in the vrp T-matrix, Eq.(|TB|). The generic term is 
of the form 

Tf, = f{Ql) [zS,-{AxC)], (HI) 

where (5+ = C + A. (Additional dependence on the rotational scalar A"^ = (7^ is a trivial 
modification of this angular decomposition.) To proceed, we expand /(Q+) in spherical 
harmonics, as in Eq.(Gl); 

fiQl) = E fM') E Y;^{Qc)Yemi^A) (H2) 

i m 

and introduce spherical components for the spin and momentum vectors, 

{Is'JS^lls,) = -V2 {Is'Jlfi, Is,) , (H3) 

t{AxC), = ^^A' {lf^',lf^"M Y,,iQA) Y.^.iQc) (H4) 

and the usual state vector expansion, 

\jm, Is {PsV) ) = J2 (^/^. ^Sz\jm) Yi,{Qa) \ls,) . (H5) 

With these substitutions one may determine the up {jl' s\Tfi\jls) matrix elements (analo- 
gous to the spinless matrix element Ti of Eq.(D7)) for the spin-orbit term (HI). The result 
involves a sum over a product of six Clebsch- Gordon coefficients, and can be written as 
the product of two Wigner (3j) symbols and two {6j} symbols, 
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{jl's\Tf,\jls) = i-iy^'GA' (2^+1) 7(2/ +1) (2/' +1) 





^1 /' e) 


1 




II 




^0 0^ 


^0 j 











(H6) 



The constraints of the (3j) and {Qj} symbols force this matrix element to be diagonal in 
1,1', and imply that the only radial components of f{Q^) in Eg. (P2|) that contribute are 
fi=i±i. Substitution of the explicit (3j) and {6j} symbols gives our final result for PsV 
matrix elements of spin-orbit (HI) type, 

This result has the overall (L ■ S ) dependence that one would expect from a spin-orbit 
force. 

We may similarly evaluate the matrix elements of the tensor terms in Eq.(16). It 
suffices to consider the two cases 

TfP = f{Ql) {[S,-AS,-A- ^-A'] + [S^.CS,-C-'^-C']) (H8) 

and 

Tjf = /(Q'_) [S,-AS,-C-'^-A-C] . (H9) 

Both tensor matrix elements have I ^ V contributions, unlike the other interactions we 
have considered. The general results in terms of Wigner (3j) and {6j} symbols are 

01/2 cl/2 2 

{3l's\TfP\jls) = A (/, + /,) v'(2/ + l)(2/' + l) 



(2, A 


1 




^0 j 







(HIO) 



and 



41 



{jl's\Tif\jls) = i-iy^'BA' (2^ + 1) V(2^ + l)(2/' + l) 





\ v A 


1 




II 




^0 0^ 


^0 0^ 











(Hll) 



Substitution for the (3j) and {6j} symbols gives the results quoted in Eqs. (23-25) in the 
text. 



42 



[1] J.Weinstein and N.Isgur, Phys. Rev. Lett. 48, 659 (1982); Phys. Rev. D27, 588 (1983); Phys. 
Rev. D41, 2236 (1990). 

[2] K.Dooley, E.S.Swanson and T.Barnes, Phys. Lett. B275, 478 (1992). 

[3] N.TSrnqvist, Z. Phys. C61, 525 (1994). 

[4] P.L.Frabetti et al. (E687 Collaboration), Phys. Lett. B351, 591 (1995). W.-S.Hou and K.- 
C.Yang, |hep-ph/9911528| v3, suggest that FSI effects may also be important in B decays. 



[5] R.Machleidt, K.Hohnde and C.Elster, Phys. Rep. 149, 1 (1987). A recent variant of this 
approach derives NN scattering amplitudes from meson exchange at the quark level; see 
D.Bartz and Fl.Stancu, |hep-ph/00060T^ . 



[6] There is an extensive literature on the derivation of the NN interaction from quark-gluon 
forces in the quark model, beginning with D.A.Liberman, Phys. Rev. D16, 1542 (1977). 
Most studies have employed resonating group methods (see for example C.S.Warke and 
R.Shanker, Phys. Rev. C21, 2643 (1980); M.Oka and K.Yazaki, Phys. Lett. 90B, 41 (1980); 
J.E.T.Ribeiro, Z.Phys.C5, 27 (1980)) and variational methods (for example K.Maltman and 
N.Isgur, Phys. Rev. D29, 952 (1984)). The Born-order quark-gluon approach described here 
was applied to NN by T.Barnes, S.Capstick, M.D.Kovarik and E.S.Swanson, Phys. Rev. 
C48, 539 (1993), and in a closely related approach by D.Hadjimichef, G.Krein, S.Szpigel 
and J.S.daVeiga, Phys. Lett. B367, 317 (1996). 

[7] M.Liischer, Comm. Math. Phys. 104, 177 (1986); 105, 153 (1986); Nucl. Phys. B354, 531 
(1991); S.R.Sharpe, R.Gupta and G.Kilcup, Nucl. Phys. B383, 309 (1992); R.Gupta, A.Patel 
and S.R.Sharpe, Phys. Rev. D48, 388 (1993); M.Fukugita et al, Phys. Rev. Lett. 71, 2387 
(1993); Phys. Rev. D52, 3003 (1995); S.Aoki et al, |hep-lat/99 110251 (Fig.5 shows the 
scattering length quoted by Aoki et al.) 



43 



[8] C.Stewart and R.Koniuk, Phys. Rev. D57, 5581 (1998); C.Michael and P.Pennanen, Phys. 
Rev. D60, 054012 (1999); P.Pennanen, C.Michael and A.M.Green, |hep-lat/990803^ . 



[9] T.Barnes and E.S.Swanson, Phys. Rev. D46, 131 (1992); see also G.-Q.Zhao, X.G.Jing 
and J-C.Su, Phys. Rev. D58, 117503-1 (1998); D.Hadjimichef, G.Krein, S.Szpigel and 
J.S.daVeiga, Phys. Lett. B367, 317 (1996); Ann. Phys. (N.Y.) 268, 105 (1998); B.Masud, 
J.Paton, A.M.Green and G.Q.Liu, Nucl. Phys. A528, 477 (1991); B.Masud, Phys. Rev. D50, 
6783 (1994). 

[10] E.S.Swanson, Ann. Phys. (N.Y.) 220, 73 (1992). 

[11] T.Barnes, E.S.Swanson and J.Weinstein, Phys. Rev. D46, 4868 (1992). 



[12] T.Barnes and E.S.Swanson, Phys. Rev. C49, 1166 (1994); ses also D.Hadjimichef, fiep- 
I ph/0006330| and |nucl-th/9912"02^ . 



[13] T.Barnes, N.Black, D.J.Dean and E.S.Swanson, Phys. Rev. C60, 045202 (1999). The Phys. 
Rev. C version has an incorrect figure (a duplicate of Fig. 4 was used as Fig. 5), and minus 
signs were deleted from the lower y-axis in both figures. The original, archived version 



iiucl-th /9902068| has the correct figures. 



[14] D.Blaschke and G.Ropke, Phys. Lett. B299, 332 (1993); K.Martins, D.Blaschke and 
E.Quack, Phys. Rev. C51, 2723 (1995); K.Martins, Prog. Part. Nucl. Phys. 36, 409 (1996); 
D.Blaschke, G.Burau and Yu.L.Kalinovsky, |nucl-th/0006071 . 



[15] C.-Y.Wong, E.S.Swanson and T.Barnes, hep-ph/9912431| (Phys. Rev. C to appear); [nucl 



th/0002034t |nucl-th/00060T2 . 



[16] G.Janssen, K.Holinde and J.Speth, Phys. Rev. C49, 2763 (1994). 

[17] R.Bockmann, C.Hanhart, O.Krehl, S.Krewald and J.Speth, Phys. Rev. C60, 055212 (1999). 
[18] E.Colton et al, Phys. Rev. D3, 2028 (1971). 
[19] N.B.Durusoy et al, Phys. Lett. 45B, 517 (1973). 

44 



[20] W.Hoogland et al, Nucl. Phys. B126, 109 (1977). 
[21] M.J.Losty et al, Nucl. Phys. B69, 185 (1974). 
[22] J.P.Prukop et al, Phys. Rev. DIO, 2055 (1974). 

[23] D.Lohse et al, Nucl. Phys. A516, 513 (1990); G.Janssen et al, Phys. Rev. D52, 2690 (1995). 
[24] F.Sassen, personal communication. 



[25] B.Ananthanarayan, G.Colangelo, J.Gasser and H.Leutwyler, hep-ph /0005297 . 
[26] S.Weinberg, Phys. Rev. Lett. 17, 616 (1966). 

[27] D.Morgan and M.R.Pennington, "Low Energy tttt Scattering." , University of Durham report 
DTP-95/28 (March 1995). 

[28] M.R.Pennington, Nucl. Phys. A623, 189c (1997). 

[29] J.F.Donoghue, C.Ramirez and G.Valencia, Phys. Rev. D38, 2195 (1988). 

[30] N.F.Mott and H.S.W.Massey, The Theory of Atomic Collisions (Clarendon Press, Oxford 
1965), Sec.VIL8.1. 

[31] T.Barnes and G.I.Ghandour, Phys. Lett. B118, 411 (1982). 

[32] S.Krewald, personal communication. 

[33] Fl.Stancu, D.M.Brink and H.Flocard, Phys. Lett. 68B, 108 (1977). 

[34] D.Mukhopadhyay and H.-J.Pirner, Nucl. Phys. A442, 605 (1985). 



[35] For 7ri(1600) see V.Dorofeev (VES Collaboration), in Proceedings of WHS99, |hep- 



ex/9905002| ; G.S.Adams et al. (E852 Collaboration), Phys. Rev. Lett. 81, 5760 (1998). For 
7ri(1400) see D.R.Thompson et al. (E852 Collaboration), Phys. Rev. Lett. 79, 1630 (1997); 
A.Abele et al. (Crystal Barrel Collaboration), Phys. Lett. B423, 175 (1998); S.U.Chung et 
al. (E852 Collaboration), Phys. Rev. D60, 092001 (1999). 



45 



[36] N.Isgur, R.Kokoski and J.Paton, Phys. Rev. Lett. 54, 869 (1985); T.Barnes, F.E.Close and 
E.S.Swanson, Phys. Rev. D52, 5242 (1995); C.Bernard et al. (MILC Collaboration), Phys. 
Rev. D56, 7039 (1997); P.Lacock et al. (UKQCD Collaboration), Phys. Lett. B401, 308 
(1997). 



[37] A.Donnachie and P.R.Page, |hep-ph/9807433|; Phys. Rev. D58, 114012 (1998) 



[38] M.L.Goldberger and K.M.Watson, Collision Theory (J.Wiley, New York 1964). 

[39] Particle Data Group, Review of Particle Properties, Eur. Phys. J. C3, 1 (1998). 

[40] The possibility of measuring this phase using Crystal Barrel bi ujtt data is being investi- 
gated by E.Klempt and B.Pick. 

[41] L.I.Schiff, Quantum Mechanics (McGraw-Hih, New York, 1968), pp.384-387. 



46 



C1 



C2 




T1 




T2 



Fig.l. The four quark interchange meson- meson scattering diagrams. 
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Fig. 2. The quark-quark T- matrix, showing three-momentum definitions. 
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Fig. 3. Theoretical contributions to the 1=2 tttt S-wave phase shift, Eq.(6), with SHO 
wavef unctions. The experimental phase shifts of Colton et al. JlSl (down triangles), Du- 
rusoy et al. 0] (up triangles, two extrapolations), Hoogland et al. (set B, diamonds), 
and Losty et al. (squares) are shown. 
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Fig. 4. Numerically evaluated 1=2 tttt S- and D-wave phase shifts with Coulomb plus linear 
plus hyperfine wavefunctions (lines), compared to experimental phase shifts (symbols as 
in Fig. 3). 
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Fig. 5. A "generalized specific heat plot" of the 1=2 tttt S-wave phase shift. The data of 
Fig.3 is shown together with a^=2 predictions: a) LGT 0], Roy Eqs. |2|, b) PCAC |2|, 
xPT [^. Meson exchange p3|j2^ and quark model (Eq.(6)) predictions are also shown. 
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Fig. 6. The low energy 1=2 tttt S-wave potential, Eq.(13). 



49 




_45 I ' ' ' ' ' ' ' ' ' 

0.5 1 1.5 2 2.5 

M.p[GeV] 



Fig. 7. The theoretical 1=2 irp S-wave phase shift with SHO 
wavef unctions, Eq.(28). 
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Fig. 8. Theoretical 1=2 np phase shifts in P-, D- and F- waves. 
^P2, ^Pi, ^Pq and J=L±1 phase shifts are shown. 
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Fig. 9. OGE spin-orbit, confining spin-orbit and OGE tensor 
contributions to the '^P2 1=2 vrp phase shift. 




Fig. 10. Spin-orbit potentials in the 1—2 np channel, Eq.(29). 
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